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We provide a theoretical estimate of the thermoelectric current and voltage over a Phenyldithiol 
molecule. We also show that the thermoelectric voltage is (1) easy to analyze, (2) insensitive to the 
detailed coupling to the contacts, (3) large enough to be measured and (4) give valuable information, 
which is not readily accessible through other experiments, on the location of the Fermi energy relative 
to the molecular levels. The location of the Fermi-energy is poorly understood and controversial 
even though it is a central factor in determining the nature of conduction (n- or p-type). We also 
note that the thermoelectric voltage measured over Guanine molecules with an STM by Poler et al. , 
indicate conduction through the HOMO level, i.e., p-type conduction. 



I. INTRODUCTION 

Electrical conduction through individual molecules 
chemically bound to gold contacts has recently been mea- 
sured using STM 0,0] j break junctions Q and nanopores 
Theoretical calculations of the current- voltage (I-V) 
characteristics have also been reported, see for example 
Ref. 0, S S 13 ■ However, a detailed quantitative com- 
parison of theoretical and experimental results is made 
difficult by two factors. Firstly, the low bias conductance 
depends strongly on the quality of the metal-molecule 
contacts Q, which is ill controlled and poorly charac- 
terized experimentally. Secondly, the conductance gap is 
determined by the location of the Fermi energy relative to 
the molecular levels - a factor that is poorly understood 
and controversial jsL Il0| - This is particularly unfortunate 
since the location of the Fermi energy is a central factor 
in determining the nature of conduction (n- or p-type) . 

The purpose of this paper is to show that a measure- 
ment of the thermoelectric voltage (see Fig. ^) can pro- 
vide new insights into electron transport and allows us 
to estimate the location of the Fermi energy relative to 
the molecular levels. Analogous to the hot point probe 
measurements commonly used to establish the p- or n- 
character of semi-conductors (see for example Ref. ^2), 
the thermoelectric voltage yield valuable information re- 
garding the location of the Fermi energy. In contrast 
to the I-V, the thermoelectric voltage is an equilibrium 
property which is easy to interpret. What makes it par- 
ticularly useful is that unlike the low-bias conductance, 
it is relatively insensitive to the quality of the contacts. 

In this paper we estimate the thermoelectric current 
and voltage over a Phenyldithiol (PDT) molecule, see 
Fig. ^ We show that the thermoelectric voltage is 
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FIG. 1: a: Proposed experimental setup to measure thermo- 
electric effects over a molecule with two contacts at different 
temperatures T\, Ti- In this paper we focus on the Phenyl- 
dithiol (PDT) molecule (1) chemically bound to both contacts 
(break junction setup) and (2) with one weak contact (STM 
measurement), b, c: Equivalent circuits defining the voltage, 
conductance and current. 



(1) easy to analyze, (2) insensitive to the contacts and 
(3) large enough to be measured. In fact, Poler et al. 
[12] measured the thermoelectric voltage (with a 20 K 
temperature difference) over a monolayer of Guanine 
molecules with an STM. In this paper we also show that 
the thermoelectric voltage provide information on where 
the Fermi energy is relative to the molecular levels, e.g., 
the sign of the thermoelectric voltage over the Guanine 
molecule indicate p-type conduction. 



II. METHOD 



'Electronic address: 
t Electronic address: 



mpaulsso@purdue.edu 



datta@purdue.edu 



The electrical transport through small conjugated 
molecules chemically bound to at least one of the con- 
tacts can be described by the Landauer formula in terms 
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of the transmission (T(E)) [l3j: 

oo 

I=^J T{E){h{E)-f 2 {E)) dE (1) 



assuming that the electronic states in the contacts are 
filled according to the Fermi distribution of the reservoirs 
(/lj fi)- If the transmission is not constant (as a function 
of energy), a difference in the Fermi distributions due to 
a temperature difference (AT = T\ — T 2 ) will drive a 
thermoelectric current at zero bias (V — 0) (also see 
footnote [lj). For molecules the transmission is often 
smooth (compared to the thermal e nerg y) which allows 
us to use the Sommerfeld expansion |14| : 
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AT (2) 



E=E, 



where T is the mean temperature of the contacts ((Ti + 
T 2 )/2). Similar expressions are available for heat trans- 
port [3. 

The thermoelectric current is usually small enough 
that we can use a linear equivalent circuit as shown in 
Fig. \I\b. The current source (see Eq. |2J) has a resistance 
that is the inverse of the low bias conductance: 



Reed et al. jTe]. In any case, our objective is to provide 
a reasonable qualitative estimate rather than an accu- 
rate quantitative value for a specific molecule. On gold 
surfaces, PDT chemisorbs forming a bond between the 
sulfur and gold. The strength of this bond is difficult to 
estimate since the precise experimental geometry is un- 
known. For an STM measurement the contact with the 
tip is also weaker than the substrate bond. We therefore 
calculate the transmission through the PDT molecule 
perfectly bound to the gold contacts (see molecule 1 in 
Fig. n a) and through a PDT molecule chemisorbed 
on only one contact with different distances between the 
molecule and contact 1 (see molecule 2 in Fig. ^a,). 
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Alternatively we could represent the effect in terms of 
a voltage source (V = I/G) with a series resistance as 
shown in Fig. ^ c: 
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As we will see below, this open circuit voltage can be 
large enough to measure and is relatively insensitive to 
the coupling to the contacts. 

To obtain estimates of the low bias conductance and 
the thermoelectric voltage we need to calculate the trans- 
mission as a function of energy for a molecule connected 
to metallic contacts. Here we calculate the transmission 
using the Non Equilibrium Green's Function (NEGF) for- 
malism using Extended Hiickel (EH) theory (see footnote 
|22|) as described in Ref. ^5|. This is a simple but suf- 
ficient approximation if we take Ef as an adjustable pa- 
rameter (see footnote |23j l to be determined from exper- 
iments since the qualitative features of the transmission, 
in our experience, is not very sensitive to the method 
(Hiickel, EH or Ab-initio) used. 



III. AN ILLUSTRATIVE EXAMPLE 

We choose to estimate the thermoelectric voltage 
across a phenyl-dithiol (PDT) molecule (Fig. [IJi) since 
it has been studied extensively after the experiment by 



FIG. 2: Logarithmic plot of the transmission as a function 
of energy. Thin solid line: Transmission calculated from the 
Extended Hiickel model through the PDT molecule strongly 
bound to both contacts. Thin dashed line: PDT weakly 
bound to one side (Au-H distance 2.9A). Thick lines: The 
transmission fitted to Lorentzian transmission peaks associ- 
ated with the molecular energy levels (see Eq. |^J. Note the 
widely different positions of the Fermi energy (Ef) suggested 
by different authors. 

Fig. |21 shows the transmission through the PDT 
molecule for perfect contacts on both sides (solid line) 
and one weak contact (dashed line). The position of 
the Fermi energy (Ef) relative to the Highest Occu- 
pied Molecular Orbital (HOMO) and Lowest Unoccupied 
Molecular Orbital (LUMO) levels is one of the main pa- 
rameters affecting the current- voltage characteristics [j| . 
As pointed out in the introduction, the position of Ef 
is difficult to estimate. The suggested positions of Ej 
range from (1) closer to the HOMO level [1 H [H Q]f , 
(2) mid gap and (3) closer to the LUMO level |[g, see 
Fig. |21 Since the thermoelectric effect depends on the 
slope of the transmission at Ef (Eq. 0J, the sign of the 
thermoelectric voltage and current will be determined by 
where Ef is located relative to the molecular levels. 

The thermoelectric voltage calculated from Eqs. 0J[Slis 
shown in Fig. |3]for a 10 K temperature difference at room 
temperature (T x = 300 K, T 2 = 310 K). Note that the 
thermoelectric voltage was calculated from a transmis- 
sion fitted to Lorentzian peaks (see Fig. |2J), this removes 
the magnification (by taking the derivative) of small nu- 
merical errors. Unless Ef is located very close to the 
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middle of the HOMO-LUMO gap the size of the thermo- 
electric voltage is of the order 0.5 — 0.1 mV. Since the 
thermoelectric voltage is proportional to AT (Eq. 0J, 
even smaller temperature differences will give a measur- 
able thermoelectric voltage. 
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FIG. 3: Thermoelectric voltage generated by a 10 K differ- 
ence in temperature (Ti = 300 K, T2 = 310 K) for different 
positions of the Fermi energy (Ef). Solid line: PDT strongly 
bound to both contacts. Dashed line: weak bond on one side. 
Thin dotted line, approximated by Eq. |S| Inset: Logarithmic 
plot of the thermoelectric current. 



The measured conductance of the PDT molecules 0] 
is smaller than the theoretical estimates by two to three 
orders of magnitude ||. This discrepancy indicate that 
the interaction strength between molecule and contacts 
is smaller than that assumed in the models (for STM 
measurements this is obviously true). It is plausible to 
believe that the molecule form a strong bond with one of 
the contacts and only interacts weakly with the second 
contact (also see footnote (24|). However, the calculated 
thermoelectric voltage with one weak contact (Fig. |3| 
dashed line) has almost the same magnitude as the per- 
fectly bound molecule (solid line) and is slightly larger 
close to the transmission peaks. 

To understand why the thermoelectric voltage is unaf- 
fected by weak contacts it is instructive to fit the trans- 
mission (Fig. [21 to Lorentzian transmission peaks: 
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t (E-etf + (T 1 +T 2 f /A 



(5) 



where e$ is the energy of the two levels and r^,^ the 
broadening by contact 1 and 2. This approximation 
is useful for the PDT molecule since the transmission 
around the HOMO-LUMO gap of PDT is mainly deter- 
mined by two levels that couple approximately equally 
to the contacts (see footnote |25|). Fitting the transmis- 
sion to Eq. fusing a least square fit (of hi(T)) gives: 
Ti = T 2 = 0.11 eV, ei = -11.86 eV and e 2 = -7.91 eV. 
For the weakly coupled PDT molecule (with a Au-H dis- 
tance of 2.9 A) the least squares fit gives (r 2 = 0.11 eV): 
Ti = 0.0042 eV, ei = -11.80 eV, e 2 = -7.85 eV. These 



fits give good approximations to the calculated transmis- 
sion (see Fig. other Au-H distances (not shown in 
figure) in the range 1 — 5 A also fit the transmission well. 

Assuming the Fermi energy to be situated in between 
the HOMO and LUMO levels and far from the levels 
(\Ef - ei )2 | >> |Ti + r 2 |) we can Taylor expand Eq. 0| 
in the energy around the midpoint of the HOMO-LUMO 
gap (E = a±a): 



87T 2 fc|T 
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E f 



ei + e 2 



AT (6) 



This shows that the thermoelectric voltage is, to first 
order, independent of the strength of the interaction with 
the contacts, see thin dotted line in Fig. [3] This result 
is not limited to the PDT molecule but valid for most 
conjugated molecules (see footnote prjjl. 

To experimentally measure the thermoelectric voltage, 
the thermoelectric current has to be larger than any leak- 
age currents. The thermoelectric current (Eq. [2J through 
the PDT molecule is shown in the inset of Fig. [21 Since 
the thermoelectric current is proportional to the trans- 
mission (T), the magnitude of the current is strongly 
dependent on the strength of the interactions between 
molecule and contacts. It is also important to realize 
that for very weak coupling to both contacts, i.e., in the 
Coulomb blockade regime, the simple decomposition of 
the transmission into HOMO and LUMO peaks (see Fig. 
EJ) might not be enough to describe the thermoelectric 
effect. However, even for a molecule which is relatively 
weakly connected to one contact (but chemically bound 
to the other contact) should give a thermoelectric current 
of the order of 1 — 100 pA, see inset in Fig. 



IV. SUMMARY AND DISCUSSION 

The position of the Fermi energy is one important in- 
gredient in understanding electrical transport through 
molecules. It is also interesting since shifting the Fermi 
energy by, (1) doping, (2) addition of side-groups, (3) 
gate field or (4) use of different contact material, can be 
used to optimize the molecular properties. Phenomena 
of this type can be probed and evaluated through a mea- 
surement of the thermoelectric voltage (coupled with I-V 
measurements) . 

In this study we have only considered a single molecule 
between the two contacts. However, unless the local 
structure of each molecule is different, the thermoelectric 
effect should only be affected trivially. It is possible that 
the detailed structure of the gold atoms bonding to the 
molecule affect the charge transfer to the molecule and 
thus the position of the Fermi energy. If this is the case, 
the variations could be mapped by STM measurements 
on different molecules on a single substrate. Another sim- 
plification in our work is that we only focus on the PDT 
molecule. Due to its simple structure the transmission 
is almost symmetric around the middle of the HOMO- 
LUMO gap. For other molecules this need not be true. 



4 



In this case, the interpretation of the thermoelectric volt- 
age still provides an estimate of the Fermi energy even 
though it is less straight forward. 

To our knowledge the only reported measurement of 
the thermoelectric voltage over a molecule was performed 
by Poler et al. 12]. Using an STM tip they measured 
the thermoelectric voltage over Guanine molecules on a 
graphite substrate. A detailed analysis of the measure- 
ment is not possible since the I-V characteristic was not 
measured in this experiment. However, the sign of the 
thermoelectric voltage indicates that the electrical trans- 
port is conducted through the HOMO level and the mea- 
sured thermoelectric voltage was similar in magnitude to 
the value estimated here for PDT ( 0.5 ± 0.01 mV at 
AT = 26.5 K). 

An intriguing possibility is to use the thermoelectric 



effect over a molecular SAM as a thermoelectric ele- 
ment. Applications of these elements include thermo- 
electric coolers (Peltier effect) and power generators. We 
have performed preliminary calculations on the efficiency 
of such elements. However, the efficiency is limited (and 
largely determined) by the lattice thermal conductivities 
of the SAM. We therefore believe that until there are 
reliable estimates of the thermal conductivity any specu- 
lations on the efficiency of such elements are premature. 

In summary, we have shown that a measurement of 
the thermoelectric voltage over a molecule is (1) simple 
to analyze, (2) insensitive to the contacts, (3) feasible and 
(4) should give valuable information, which is not readily 
accessible through other experiments, on the location of 
the Fermi energy relative to the molecular levels. 
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